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ABSTRACT 

We present FUSE spectroscopy of SBS 0335-052, the second most metal-deficient 
blue compact dwarf (BCD) galaxy known (log O/H = -4.70). In addition to the Hi 
Lyman series, we detect Cll, Nl, Nil, Oi, Sin, An and Fen absorption lines, mainly 
arising from the extended Hi envelope in which SBS 0335-052 is embedded. No H2 
absorption lines are seen. The absence of diffuse H2 implies that the warm H2 detected 
through infrared emission must be very clumpy and associated with the star-forming 
regions. The clumps should be denser than ~ 1000 cm -3 and hotter than ~ 1000 K and 
account for > 5% of the total Hi mass. Although SBS 0335-052 is a probable young 
galaxy, its neutral gas is not pristine. The metallicity of its neutral gas is similar to that 
of its ionized gas and is equal log O/H ~ -5. This metallicity is comparable to those 
found in the Hi envelopes of four other BCDs with ionized gas metallicities spanning 
the wide range from log O/H = -4.8 to log O/H = -3.8, and in Lya absorbers, fueling 
the speculation that there may have been previous enrichment of the primordial neutral 
gas to a common metallicity level of log O/H ~ -5, possibly by Population III stars. 



Based on observations obtained with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer. FUSE is 
operated for NASA by the Johns Hopkins University under NASA contract NAS 5-32985. 
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Subject headings: galaxies: irregular — galaxies: abundances — galaxies: ISM - 
galaxies: individual (SBS 0335-052) 

1. Introduction 

The Blue Compact Dwarf (BCD) galaxy SBS 0335-052 (Izotov et al. 1990), with an ionized 
gas oxygen abundance of 12 + log O/H = 7.30 (Melnick et al. 1992; Izotov et al. 1997, 1999) 
[Zq/23, if the recent lower determination of the Oxygen abundance of the Sun 12 + log O/H = 
8.66 by Asplund et al. (2004) is adopted], is one of the most metal-deficient star-forming galaxies 
known, just after its dwarf irregular companion SBS 0335-052W (12 + log O/H = 7.22 or Z Q /28, 
Lipovetsky et al. 1999) and the BCD I Zw 18 (12 + log O/H = 7.17 or Z /31, Izotov et al. 
1999). Because of its extremely low metallicity and its high neutral hydrogen content (Pustilnik 
et al. 2001), SBS 0335-052 is an excellent candidate for being a young galaxy, just like I Zw 18. 
The latter has been shown by Izotov & Thuan (2004) to be a bona fide young galaxy. Resolving 
I Zw 18 into stars with deep Hubble Space Telescope (HST) /Advanced Camera for Survey images 
and constructing its color-magnitude diagram (CMD), Izotov & Thuan (2004) have shown that 
the galaxy does not contain red giant stars and that the most evolved stars in it are not older than 
500 Myr. Because of its larger distance (54.3 Mpc, Thuan et al. 1997, hereafter TIL97, instead 
of 15 Mpc for I Zw 18), SBS 0335-052 cannot be resolved into stars by HST, so its age cannot be 
determined directly by CMD analysis. However, using optical colors of the low surface brightness 
component together with evolutionary synthesis models, TIL97 and Papaderos et al. (1998) have 
constrained the age of the underlying stellar population to be less than ~ 100 Myr. Near-infrared 
(NIR) colors after correction for gaseous emission are consistent with a stellar population not older 
than 4 Myr, with a possible contribution to the NIR light from an evolved stellar population not 
exceeding ~ 15% (Vanzi et al. 2000). Ostlin &: Kunth (2001) have derived from surface photometry 
considerably larger ages for SBS 0335-052. However these large ages may be the consequence of 
their not subtracting out the large ionized gas contamination from the surface brightness profiles 
used. Recently, Pustilnik et al. (2004) have also found a young age for SBS 0335-052. Depending 
on the adopted star formation history of the BCD, they derive an age between 100 Myr and 400 
Myr. In any case, SBS 0335-052 constitutes an excellent nearby laboratory for studying massive star 
formation and its interaction with the ambient interstellar medium (ISM) in a very metal-deficient 
environment. Most of the star formation in SBS 0335-052 occurs in six super-star clusters, roughly 
aligned in the southeast-northwest direction, with ages < 25 Myr, within a region of ~ 2" or 520 
pc in size (TIL97). 

VLA observations have shown the presence of a massive (2.2 xlO 9 M ) and extended (66 
kpcx22 kpc) neutral hydrogen gas envelope around SBS 0335-052 (Pustilnik et al. 2001). This 
Hi envelope absorbs all the Lya photons emitted by the young massive stars. Thuan & Izo- 
tov (1997) have obtained a HST/Goddard High Resolution Spectrograph (GHRS) spectrum of 
SBS 0335-052 which reveals a broad damped Lya absorption with a high Hi column density 
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iV(Hi) = (7.0±0.5)xl0 21 cm' 2 , the largest known in a BCD. The fTST/GHRS spectrum also 
shows absorption lines of several heavy elements such as O I, Si n and S II. Assuming that these 
absorption lines are not saturated, Thuan & Izotov (1997) found that the derived heavy element 
abundance in the neutral gas to be considerably lower than that of the ionized gas. Furthermore, 
Vanzi et al. (2000) have detected with infrared spectrophotometry molecular hydrogen in emission 
in SBS 0335-052. The presence of a neutral gas envelope around SBS 0335-052 and the detection 
of molecular hydrogen in emission make this galaxy an ideal target for FUSE spectroscopy. Using 
SBS 0335-052 as a source of UV light shining through the H I envelope, a FUSE absorption spec- 
trum will allow us to determine independently the heavy element abundances in the neutral gas of 
the galaxy and check for the presence of H2 in it. 

We describe the FUSE observations in Section 2. In Section 3, we set upper limits on the 
amount of diffuse H2. In Section 4, we show that the warm H2 detected through infrared emission 
must be very clumpy. We derive the column densities of the interstellar ionic species in the H 1 
gas, and compare the heavy element abundances in the neutral and ionized gas in Section 5. We 
summarize our findings in Section 6. 

2. Observations and data analysis 

SBS 0335-052 has been observed with FUSE (Moos et al. 2000) through the 30"x30" LWRS 
large entrance aperture on 2001, September 26 for a total integration time of 24.4 ksec. The data 
have been processed with version 1.8.7 of the CALFUSE pipeline. The eight separate exposures have 
been aligned and coadded. This results in a set of four independent spectra, two spectra for the 
two long wavelength LiF channels covering the wavelength range ~ 1000-1187 A, and two spectra 
for the two short wavelength SiC channels, covering the wavelength range ~ 900-1100 A. After 
co-addition, the S/N ratio per resolution element and per channel is estimated to be about 2 below 
1000 A and about 4 above 1000 A. 

To derive the spectrum, special attention was paid to the background residual which is not 
negligible for this faint target. The background has been estimated both on spectra obtained off- 
target through other apertures and at the bottom of wide saturated lines. The background level is 
found to be between and 3 x 10~ 15 erg cm~ 2 s _1 A -1 . The line spread function has been estimated 
with the unresolved Galactic H2 absorption lines. The resolution is found to be R = AA/A ps 12000. 

The resulting spectrum of SBS 0335-052 shifted to its rest-frame is shown in Fig. 1. It 
exhibits two absorption-line systems at two different radial velocities. The first system at nearly 
zero radial velocity is attributed to the interstellar clouds in the Milky Way. The second system 
at ~4050 km s" 1 is attributed to the ISM in SBS 0335-052, as the Hi velocity of the BCD is 
4057±5 km s _1 (Pustilnik et al. 2001), in good agreement with the optical velocity obtained 
by Izotov et al. (1997). Absorption lines from the atomic and ionic species Ni, Pn, Fen, and 
Nil! are observed in the Milky Way, while for SBS 0335-052, in addition to the H 1 Lyman series, 
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absorption lines from the atoms and ions Cn, Ni, Nil, Oi, Sin, An, and Fen are seen. In 
Fig. 1, the lines arising from SBS 0335-052 are labeled while those belonging to the Milky Way 
are indicated by tickmarks below the FUSE spectrum. For the characteristics of the electronic 
transitions of H2, we have used the wavelengths and oscillator strengths tabulated by Abgrall et 
al. (1993a,b) respectively for the Lyman and the Werner systems, and the inverses of the total 
radiative lifetimes tabulated by Abgrall et al. (2000). For atomic lines, we have used the data for 
resonance absorption lines tabulated by Morton (2003). 



3. Upper limits on the diffuse H 2 content of SBS 0335 052 

Molecular hydrogen is known to be present in SBS 0335-052 as Vanzi et al. (2000) have 
detected several H2 emission lines in their low-resolution near-infrared spectrum of the star-forming 
region. These lines observed in the K band are generally consistent, within the errors, with both 
thermal and fluorescent excitation by the strong UV field. It is thus of great interest to check 
whether there are H 2 absorption lines in the FUSE spectrum of SBS 0335-052. We found no H 2 
line at the radial velocity of the BCD (Fig. 1). We obtained the following 3 a upper limits for the 
H 2 column density: 1.2xl0 17 cm" 2 , 1.4xl0 17 cm" 2 , 7.9xl0 16 cm" 2 , 3.4xl0 17 cm" 2 and 1.5xl0 17 
cm" 2 respectively for the J= 0, 1, 2, 3 and 4 levels. To set an upper limit on the total H 2 column 
density, we follow the procedure of Aloisi et al. (2003). We add the two upper limits for the 
J=0 and J=l levels. For the higher J levels, we assume a temperature (T = 500 K and T = 1000 
K), which determines the level populations, and then normalize the level populations so as to be 
consistent with the measured upper limits. We obtain upper limits for the total H 2 column density 
of 3.7xl0 17 cm" 2 for T = 500 K and 5.5 xlO 17 cm" 2 for T = 1000 K. These upper limits for H 2 
are more than two orders of magnitude larger than those established for the BCDs I Zw 18 (Aloisi 
et al. 2003; Lecavelier des Etangs et al. 2004) and Mrk 59 (Thuan et al. 2002) because of 
the larger distance of SBS 0335-052, its fainter apparent magnitude and its lower signal-to-noise 
ratio FUSE spectrum. They correspond to column densities above which the damping wings start 
to broaden appreciably the H 2 absorption lines. With a H 1 column density of 7.0 xlO 21 cm" 2 as 
determined by the JfST/GHRS damped Lya profile (Thuan & Izotov 1997), this corresponds to 
a fraction of H 2 molecules /(H 2 ) = 1.1 xlO" 4 (T = 500 K) and /(H 2 ) = 1.6xl0" 4 (T = 1000 K) 
where the fraction is defined as /(H 2 ) = 2n(H 2 )/[2n(H 2 )+n(H 1)], and where n(H 2 ) and n(H 1) are 
the number densities of H 2 and H 1. 

Are the H 2 upper limits established by FUSE consistent with the expected fraction of H 2 
molecules in SBS 0335-052? This fraction can be estimated in the following manner. As discussed 
by Vidal-Madjar et al. (2000), at equilibrium, when the formation of H 2 on dust equals its 
destruction rate by UV photons, then: 

1 ' r> \ -1 rA/Yirm 

(1) 



7V(HI) 
cm -2 



Here F^ is the flux at lOOOA of UV H 2 dissociating radiation at radius Rq from the central ionizing 
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clusters in SBS 0335-052, Z is the metallicity equal to Z & /23, and iV(H i) is the H I column density. 
From the H I map of Pustilnik et al. (2001), we estimate i?o = 10", or 2.5 kpc for a H I column 
density level of 7.2xl0 20 cm™ 2 . This corresponds to a H I number density n(H i) ~ 0.1 cm -3 . Since 
the density of the H I is increasing inwards, we adopt for our calculations a mean value ro(H i) = 1 
cm™ 3 . 

The observed UV flux at lOOOA is 1.7xl0" 14 erg s™ 1 cm™ 2 A™ 1 (Thuan & Izotov 1997). 
Correcting for interstellar extinction using A(1000A) = 5.742A y (Cardelli et al. 1989) and A v = 
0.155 mag (Schlegel et al. 1998), we obtain F cor = 3.8 xlO™ 14 erg s™ 1 cm™ 2 A™ 1 . Then 

F ^= F cor(§A = 1.7 x KT 5 erg s" 1 cm™ 2 A" 1 , (2) 

where D = 54.3 Mpc is the distance of SBS 0335-052. Finally, from Eq. 1 we derive / (H2) = 
l.lxlO™ 10 , more than six orders of magnitude below and fully consistent with the upper limits 
derived from the FUSE spectrum. 



4. A clumpy interstellar medium 

Is such a low fraction of H2 consistent with the detection of near-infrared (NIR) H2 emission 
lines in the star- forming region of SBS 0335-052 by Vanzi et al. (2000)? We note that such a 
question has also arisen in the context of a circumstellar disk observed by Lecavelier des Etangs 
et al. (2001) where H2 is detected in the infrared, but where no UV absorption lines are seen. 
To address the question for SBS 0335-052, we have carried out calculations specifically for the 
H2(1,0)S(0) line with the wavelength of 2.223 [im. The results for the other detected H2 lines 
should not be too different. Vanzi et al. (2000) have obtained a flux of 1.4xl0™ 16 erg s™ 1 cm™ 2 
for the 2.223 //m line in SBS 0335-052. Assuming that the NIR emission is the result of fluorescent 
excitation, the luminosity of the H2(1,0)S(0) line from the region with radius Ro is 

L tot = n(m)f(R 2 )QhuV, (3) 

where hv = 9xl0 -13 erg is the energy of the photon with wavelength 2.223/mi, and V=4itRq/3 = 
1.92xl0 66 cm™ 3 is the volume of the region with radius Rq. Following Black Sz Dalgarno (1976) 
and using the UV flux Fr given by Eq. 2, we calculate a cascade entry rate Q = 1.55 xlO™ 11 
cm™ 3 s™ 1 . We obtain from Eq. 3 the luminosity of the 2.223/um line to be L tot =1.82x 10 33 erg 
s™ 1 . The spectrum of Vanzi et al. (2000) was obtained through a l"xl'/5 aperture, which corre- 
sponds to a volume Vy anz i = 1.38xl0 64 cm -3 . Thus, Vyanzi/V = 7.2xl0™ 3 and we would predict 
£predirfed=l-31xl0 31 erg s™ 1 , corresponding to a flux at Earth of 2.4xl0™ 23 erg s™ 1 cm™ 2 . This is 
~ 10 times lower than the observed flux. 

The large flux discrepancy comes very probably from our assumption that the interstellar 
medium of SBS 0335-052 is uniform. In all likelihood, the NIR H 2 emission comes, not from an 
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uniform low-density neutral medium, but from dense clumps. For such dense clumps, self-shielding 
of electronic transitions, which becomes important starting at iV(H2) ~ 10 14 cm -2 (Black & van 
Dishoeck 1987), further increases /(H2). Although the UV pumping decreases, thermal excitation 
of the H2 vibrational states becomes important since it scales as n(H)xn(H2). Following Black 
& van Dishoeck (1987), we find that we can reproduce well the observed flux of the H 2 (1,0)S(0) 
2.223 fim line if it comes from an ensemble of clumps with a total mass ~ 10 8 M , or about 5% 
of the total H I mass in SBS 0335-052 of 2.1 xlO 9 M Q (Pustilnik et al. 2001), and if each clump 
has a number density n(H 1) = 5xl0 3 cm -3 and a temperature T = 10 3 K. For such a clump to 
be stable against gravitational collapse, its mass has to be less than its Jeans mass which is ~ 
5xl0 4 Mq. This means that there should be at least ~ 2000 of these clumps. The fraction of H2 
molecules in such clumps is / (H 2 ) = 5xl0 3 x 6.8 x 10~ n = 3.4 x 10~ 7 . The total number of H 2 
molecules in such clumps is 7V tot = /(H^M/mn = 4.1xl0 58 , where mn = 1.673xl0~ 24 g is the 
mass of the hydrogen atom. Then the luminosity emitted in the H 2 (1,0)S(0) line is L=Nt tn(H 
i)g/u/=4.5xl0 37 erg s _1 , where q ~ 10" 10 exp(-6000/T) cm 3 s 1 is the rate of collisional excitation. 
The flux at Earth is then 1.4xl0~ 16 erg s _1 cm~ 2 in very good agreement with the observed flux 
for the H2(1,0)S(0) line at 2.223 /mi. Because of the exponential temperature dependence of the 
rate of collisional excitation, the above estimate is very sensitive to the adopted temperature for the 
gaseous clumps. This temperature cannot be lower than ~ 1000 K, since for lower temperatures, 
the total mass of the clumps needed to account for the observed H2(1,0)S(0) 2.223 /mi flux would 
exceed the H 1 mass of SBS 0335-052. 

Our FUSE observations are not sensitive to such a clumpy H2 distribution. They can only 
probe diffuse H2 along the line of sights to the several thousands of UV-bright massive stars in 
SBS 0335-052. As in the cases of the BCDs I Zw 18 (Aloisi et al. 2003; Lecavelier des Etangs 
et al. 2004) and Mrk 59 (Thuan et al. 2002), the absence of diffuse H 2 in SBS 0335-052 can 
be explained by the combined effects of a low H 1 density, a scarcity of dust grains on which H2 
molecules can form, and a large UV flux which destroys molecules. 

5. Heavy element abundances 

5.1. Column densities 

We now derive the column densities of the heavy elements in the neutral H 1 envelope of 
SBS 0335-052 by fitting the profiles of the metal absorption lines. Column densities have been 
calculated by profile fitting using the Owens procedure developed by Martin Lemoine and the FUSE 
French team. This code returns the most likely values of many free parameters like the Doppler 
widths and column densities through a x 2 minimization of the difference between the observed and 
computed profiles. The latest version of this code is particularly suited to the characteristics of 
FUSE spectra. For example, it allows for a variation of the background level, for an adjustable line 
spread function as a function of the wavelength domain, and for shifts in wavelength scale. These 
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are taken as free parameters which depend on the wavelength region and are determined by a x 2 
minimization. 

The derived column densities with their error bars are given in Table 1. We will be quoting 
1 a error bars throughout the paper. These have been estimated as usual by considering the A% 2 
increase of the \ 2 °f the fit (see Hebrard G., et al. 2002, for a full discussion of the fitting method 
and error estimation with the Owens code). These error bars include the uncertainties in the 
continuum, the intrinsic line widths (6), the background residual and the instrumental line spread 
function. It is important to remark that the use of different lines with different oscillator strengths 
of the same species allows us to constrain reasonably well all these quantities. In particular, the 
line width is well constrained by strongly saturated lines. On the other hand, the instrumental line 
spread function is mainly constrained by the very narrow lines from H2 in the Milky Way. The 
final results are obtained from a simultaneous self-consistent fit to all the data. In this final fit, 
the continuum, the background residual, the instrumental line spread function, and the physical 
parameters of the absorption lines are free parameters. They are estimated by determining the 
best x 2 m the parameters' space. 

For each wavelength domain, we use simultaneously the data from two different channels: the 
SiC 1 and SiC 2 channels for wavelengths between 940 and 990 A, and the LiF 1 and LiF 2 channels 
for wavelengths above ~ 990 A (see Sahnow et al. 2000, for a description of the FUSE channels). 
Because of the low signal-to-noise ratio of the SBS 0335-052 FUSE spectrum, only lines that are 
more or less saturated are clearly detected. This explains the relatively large error bars on the 
estimated column densities. However, for most species detected in SBS 0335-052, we have for each 
species several absorption lines from different transitions with different oscillator strengths. This 
allows us to obtain reliable column densities despite the saturation of the lines. For example, several 
1 lines are detected. The 1 column density is mainly constrained by the profiles of four strong 
transitions with rest wavelengths A = 929.5 A, 950.9 A, 988.8 A, and 1039.2 A. For Ni, we detect 
the well-known triplet at 1034-1035 A. The ions Nil and Nil* have transitions in the wavelength 
region around 1084-1085 A which is usually unusable because of its location in the hole between the 
two LiF detectors (see e.g., Lecavelier des Etangs et al. 2004, in the case of IZwl8). However, 
here the large redshift of SBS 0335-052 moves that region out of the hole and allows clear detections 
of these transitions in the LiF lb and LiF 2a channels. Sin is detected only in the Ao = 1020.7 A 
transition, the line at Ao = 989.9 A being blended with a strong Niii line. However with an optical 
depth of ~2, this line does not suffer from much saturation and allows a reliable column density 
determination, albeit with large error bars. On the other hand, C 11 which is also detected only in 
one transition at Ao = 1036.3 A, is strongly saturated. This makes its column density determination 
very sensitive to systematic errors and can explain its deviation from the abundance pattern found 
for other species (see Section 5.2). 

With the exception of Hi for which the strong damping wings in the Lyman (3 line allow a 
very accurate determination of the column density, Fe 11 is the species for which we have obtained 
the highest precision (Table 1). Fen has indeed nine clearly detected lines at Aq =1055 A, 1063 A, 
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1064 A, 1082 A, 1097 A, 1122 A, 1125 A, 1143 A and 1145 A. Those lines have optical depths ranging 
from ~2 to ~30, which allows a good determination of the Fe n column density by co-adding the 
nine separate profiles to improve the signal-to-noise ratio. Fig. 2 shows the fit to the co-added 
continuum and Fen line profile. 

To better constrain the column densities, we have also used, in addition to the FUSE data, 
the low resolution HST spectra obtained with the GHRS spectrograph and the G160M and G140L 
gratings by Thuan & Izotov (1997). These spectra show absorption lines of Cn (1334.5 A), Oi 
(1302.2 A) and Sin (1190.4, 1193.3, 1260.4 and 1304.4 A). 

For the species which are detected in both FUSE and HST spectra (Oi, Cn and Sin), the 
column densities given in Table 1 are derived from a simultaneous fit of both spectra. This assumes 
that both instruments cover the same spatial regions in the BCD, although the observations have 
been obtained through very different aperture sizes, 30" in the case of FUSE, and 2" for HST. 
This assumption is reasonable as the brightest star-forming clusters in SBS 0335-052 are confined 
within a very compact region (< 2", TIL97), and this region is fully sampled by both FUSE and 
HST apertures. That this is indeed the case is confirmed by the fact that the H I column densities 
derived independently from the FUSE and HST spectra are very similar: the FUSE column density 
is log iV(Hl) = 21.86±g;g|, while the HST column density is log iV(Hi) = 21.85 ± 0.05. The heavy 
element column densities derived from the FUSE spectra alone and from the FUSE+HST spectra 
are consistent with each other within the errors. But by using the FUSE and HST spectra together, 
we are able to decrease the error bars by a factor of about 2. For example, the derived column 
density of C n using only the FUSE spectrum is log N(C n) (FUSE) = 17.68t?;f as compared to 
log N(C n) (HST+FUSE) = 17.59lg;^. 

The simultaneous fit of a large set of lines with different oscillator strengths for the same 
species, and of many different species with the same intrinsic line width allows us to constrain 
the line broadening parameter (or Doppler width) b. In other words, a single b value is derived 
which gives the best simultaneous fit to the profile of every detected line of each species. As an 
example, Fig. 3 shows the fit to the O I A 1039 line. From the FUSE spectrum alone, we obtain 
b = U.OtH km s_1 - Thc 

error bar on the b value includes the uncertainties in the continuum 
level, the background residual, the instrumental line spread function and the column densities of 
the different species detected at the redshift of SBS 0335-052. If the FUSE and HST spectra are 
taken together, we obtain b = 11.3*\'2 km s _1 . The two b values are in good agreement and thus 
b is well constrained. 

A b value of 12 km s" 1 would correspond to a velocity dispersion a of only b/ y/2 = 8.5 km s -1 , 
somewhat smaller than the value of ~ 15 km s _1 derived from the Hi VLA map of SBS 0335-052 
(Pustilnik et al. 2001), although the two values are not directly comparable: the radio observations 
measure the velocity dispersion within a region equal to the VLA beam size of 20'.'5x 15'/0, while 
the FUSE observations probe only the star-forming region of ~ 2" in size. Note that if b is 
underestimated, the column densities derived from very saturated lines would be overestimated. 
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For instance, we estimate that if the real value of the b parameter is 16kms , the column density 
of Oi would be overestimated by ~1 dex. As for Cn, the species for which the column density 
estimate is most sensitive to the b value, it would be overestimated by ~1.3 dex. Other species 
are less sensitive to the exact value of b, making the determination of their column densities more 
robust. 

5.2. Modeling 

We use the CLOUDY code (Ferland 1996; Ferland et al. 1998, version c90.05) to construct 
a photoionized H n region model that best reproduces the optical nebular emission-line intensities 
observed in SBS 0335-052 (Izotov et al. 1997). By comparing the H n column densities predicted 
by CLOUDY with the (H I + H n) column densities observed by FUSE, we will be able to deduce 
the relative amounts of heavy elements in the H I and H n gas. 

We consider a spherically symmetric ionization-bounded H II region model. The calculations 
are stopped in the zone away from the ionizing stars where the temperature drops to 2000 K. The 
ionization in this zone is very low and it is taken to be the outer edge of the H II region. Several 
input parameters need to be set. For a distance of 54.3 Mpc and using the aperture-corrected H/3 
flux from Izotov et al. (1997), we set the H/3 luminosity to be L(H/3) = 7.9 x 10 40 erg s" 1 , and 
the number of ionizing photons to be iV(Lyc) = 1.7 x 10 53 s -1 . We use the Kurucz (1991) stellar 
atmosphere models and adopt the effective temperature of the ionizing stellar radiation to be T e g 
= 50000 K, a typical value for low-metallicity high-excitation H n regions. For the inner radius of 
the H II region, we adopt R m = 1.0 x 10 19 cm. The chemical composition of the H n region is set 
by the observed element abundances derived from optical spectroscopy of SBS 0335-052 (Izotov et 
al. 1997), except for the carbon, silicon and phosphorus abundances. For carbon we adopt log 
C/O = -0.83, and for silicon log Si/O = -1.60 (Izotov & Thuan 1999). As for phosphorus, we 
adopt the solar log P/O = -3.42 of Grevesse & Noels (1996). The adopted abundances are shown 
in Table 2. 

We run the CLOUDY code varying the filling factor / and the electron number density N e in 
order to obtain the best agreement between the predicted and observed [O n] and [O m] emission 
lines. We find that if N e varies in the ranges 7-13 cm~ 3 and / in the range 0.17 - 0.29, then we have 
adequate agreement between the observed and model line intensities. The best model is found for 
a filling factor / = 0.22 and an electron number density N c = 10 cm~ 3 . The optical line intensities 
of the best model are shown in Table 3, where we compare them with the observed ones (Izotov 
et al. 1997). The error bars of the model line intensities (and those of the model parameters in 
Tables 1 and 4) correspond to the maximum range of their values. There is good general agreement 
between observed and model line intensities, giving us confidence that the photoionization model 
is correct. The predicted total hydrogen (neutral and ionized) column density in the H II region is 
N(H I + H n) = 4.93 x 10 21 cm" 2 or log iV(H I + H n)= 21.69. As expected for a H n region 
model, most of the hydrogen is ionized: iV(H n) = 4.89 x 10 21 cm~ 2 . The column density of the 
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neutral hydrogen is two orders of magnitude lower. The model also predicts a column density of 
molecular hydrogen made via H~ N(H.2) = 1.91 x 10 11 cm" 2 , more than six orders of magnitude 
lower than our observational upper limit. 

We now compare the column densities of heavy elements predicted by CLOUDY with those 
derived from the FUSE and HST spectra. The CLOUDY predicted column densities are given in 
Table 4. Those relevant to the FUSE observations are also repeated in column 3 of Table 1. In 
Table 4, x is the radially averaged ratio of the ion number to the total number of a particular 
element, e.g., N (Fe + ) / N (Fe) . The predicted log N(X) of species X is then derived as log N(X) 
= log N(U i+H n) + log X+7H - log x(X +i ). The column densities listed in Table 1 are plotted 
with their 1 a error bars as filled circles in Fig. 4. For comparison, the corresponding CLOUDY 
predicted column densities are plotted as open circles. There are two features that should be noted. 
First, the observed and calculated relative variations from ionic element to ionic element show the 
same pattern. While this simply reflects the relative cosmic abundances of the considered species, 
retrieving that pattern from the data does show the consistency of the derived column densities for 
different species despite the low signal-to-noise ratio of our FUSE spectra. 

Second, the majority of the calculated H n column densities are systematically smaller by 
more than one order of magnitude than the observed column densities. This implies that the 
FUSE column densities arise mainly in the H I neutral envelope, and not in the gas ionized by 
young stars. 

Fig. 5 shows the ratio of the column densities of all ions relative to that of O I. As for the 
column density of O I, it is given relative to that of H I. There is general good agreement between 
the ratios of the FUSE column densities (filled circles) and the ratios of the elements abundances 
in the H n region (stars). The only exception is the N(C n)/iV(0 i) ratio which is a factor of 30 
larger than the C/O ratio in the H n region. As discussed in Section 5.1, the C n column density 
is derived from a single strongly saturated line and is probably overestimated. 

Of special interest is the N(0 i)/N(R i) ratio. The ionization potential of O I is very similar to 
that of Hi, and there is an efficient charge exchange between On and Hi. Hence there is a strong 
coupling between the oxygen and hydrogen ionization fractions. In the neutral gas, the Oi/Hi 
ratio can therefore be considered as a very good proxy for the O/H ratio and hence it is a very 
good tracer of metallicity. We obtain log N(0 i)/iV(H i) = -5.04±0.55 or [O i/H i] = -1.70. This 
is to be compared with log O/H = -4.70±0.01 for the ionized gas in SBS 0335-052 (Izotov et al. 
1999). Thus, within the errors, the metallicity of the neutral gas in SBS 0335-052 is comparable 
to that of its ionized gas. The much lower metallicity of the neutral gas of SBS 0335-052 obtained 
by Thuan & Izotov (1997) from HST observations is not correct as it is based on the erroneous 
assumption that the O I, Si n and S n lines are not saturated. 

How does the metallicity of the H I gas in SBS 0335-052 compare with those determined in 
the neutral ISM of other BCDs? Such a determination has been made for four other BCDs. For 
the most metal-deficient BCD known, I Zw 18 with a log (0/H)j of the ionized gas equal to - 
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4.82±0.01 (Izotov & Thuan 1999), Aloisi et al. (2003) have obtained log (0/H) n of the neutral 
gas = -5.37±0.28 while Lecavelier des Etangs et al. (2004) derived the higher value log (0/H) n = 
-4.7±0.35. The discrepancy between the two values may be due to the fact that, to derive O/H, 
Aloisi et al. (2003) use the strongly saturated O I line at Ao = 1039 A which, at the redshift of I 
Zw 18, is blended with terrestrial airglow. For I Zw 36 = Mrk 209 with log (O/H); = -4.23±0.01 
(Izotov & Thuan 1999), Lebouteiller et al. (2004) derived log (0/H) n = -4.5+^. For the BCD 
Mrk 59 with log (O/H); = -4.01±0.01 (Izotov & Thuan 1999), Thuan et al. (2002) have derived 
log (0/H)„ = -5.0±0.3. As for the BCD NGC 1705 with log (O/H); = -3.79±0.05 (Lee & Skillman 
2004), Heckman et al. (2001) have derived log (0/H) n = -4.6±0.3. Thus if we disregard the low 
value of Aloisi et al. (2003), the ionized gas metallicity of the five different BCDs observed by 
FUSE span a log (O/H); range from -4.82 to -3.79, i.e. a 1.03 dex spread. On the other hand, the 
metal content in their H I envelopes varies in the range log (0/H) n = -5.0 - -4.5, or [0/H] n = -1.7 
- -1.2, using the revised solar abundance (O/H) = -3.34 of Asplund et al. (2004), i.e. a 0.5 dex 
spread. 

At first glance, the metallicity spread of the neutral gas component appears to be only slightly 
lower than that of the ionized gas component. However a real difference appears once the error bars 
of the metallicity measurements are considered. The metallicities of the ionized gas, log (O/H);, 
are determined from emission-line spectra with a precision of 0.01 dex (0.05 dex for NGC 1705), 
while those of the neutral gas, log (0/H) n , derived from absorption spectra, have considerably 
larger error bars varying from 0.3 dex for Mrk 59 and NGC 1705 to 1.0 dex for I Zw 36. These 
large error bars are in fact responsible for the apparent large metallicity spread of the neutral gas. 
If we take into account the error bars, the log (O/H); values do vary over a 1.0 dex range while 
the log (0/H) n measurements are consistent with a single constant value. The spread of the log 
(0/H) n measurements around their mean value < log(0/H) n >= —4.8 can be evaluated by the 
sum of their squared normalized differences which, for five measurements, follow a x\ l aw with four 
degrees of freedom. We obtain a very low dispersion ^ [(log(0/H) n — < log(0/H) n >)/a] 2 = 1.40. 
There is more than 84% chance that \\ is larger than this value, supporting our contention that 
all five measurements are consistent with a single value of log (0/H) n . 

The fact that the most metal-deficient BCD (I Zw 18) has about the same log (0/H) n as 
the most metal-rich BCD (NGC 1705) studied thus far by FUSE suggests that the metallicities of 
the neutral and ionized gas are unrelated. As discussed in the Introduction, the two most metal- 
deficient BCDs known, I Zw 18 and SBS 0335-052 have been suggested to be bona fide young 
galaxies, with an age not exceeding 500 Myr [see Izotov & Thuan (2004) for a discussion of the 
age of I Zw 18 and TIL97, Papaderos et al. (1998) and Pustilnik et al. (2004) for a discussion of 
the age of SBS 0335-052]. Thus the H I envelopes that surround these two BCDs may be expected 
to be truly primordial, i.e. not to contain heavy elements. Our study here of SBS 0335-052 and 
previous work on I Zw 18 show that this is clearly not the case, that the H I envelopes in these 
two BCDs have been enriched previously to a level [0/H] n = -1.7 - -1.4, perhaps by Population 
III stars. This speculation appears to be supported by the oxygen abundances derived by Telfer 
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et al. (2002) for the intergalactic medium using ultraviolet absorption lines in Lya absorbers. 
Those authors derive a [O/H] range of -2.0 to —1.1 [we have corrected their values using the O 
solar abundance of Asplund et al. (2004)], which overlaps well with that derived thus far for the 
neutral ISM in five BCDs. 



6. Conclusions 

We present FUSE far-UV spectra of the second most metal-deficient Blue Compact Dwarf 
(BCD) galaxy known, SBS 0335-052 with an ionized gas oxygen abundance log O/H = -4.70±0.01. 
Because this galaxy is one of the best candidates for being a young galaxy in the local universe, with 
an age not exceeding ~ 100 Myr, the study of the abundances in its neutral hydrogen envelope is 
of special interest as it can shed light on the possible metal enrichment of the intergalactic medium 
by Population III stars several hundred million years after the Big Bang. We have obtained the 
following results: 

1. No H2 absorption lines are detected. We set a 3a upper limit for diffuse H2 in SBS 0335-052 
of 3.7xl0 17 cm~ 2 for T = 500 K and 5.5xl0 17 cm" 2 for T = 1000 K. The absence of diffuse H 2 is 
due to the combined effects of a low H 1 density, a large UV flux which destroys H2 molecules and 
a low metallicity which makes grains on which to form H2 molecules scarce. 

2. The non-detection of diffuse H2 contrasts with the detection of H2 emission lines in near- 
infrared spectra of SBS 0335-052. This implies that the detected H2 must be clumpy, and associated 
only with the dense star-forming regions. These clumps should be denser than ~ 1000 cm" 3 and 
have a temperature greater than ~ 1000 K. For such dense clumps, H2 self-shielding can protect 
the H2 molecules from dissociating UV-photons. Each clump has a mass less than ~ 5xl0 4 Mq. 
There should be at least ~ 2000 clumps, accounting for > 5% of the total H 1 mass. 

3. The O1/H1 ratio in the H 1 gas (log 7V(0 i)/N (H 1) = -5.04±0.55 is similar, within the 
errors, to the O/H ratio derived from optical emission lines in the Hn region of SBS 0335-052 (log 
O/H = -4.70±0.01), i.e. the metallicity of the neutral gas is comparable to that of the ionized gas. 
Thus although SBS 0335-052 is a probable young galaxy, its neutral envelope is not pristine. 

4. There exists a noteworthy similarity between the O/H ratio observed in the H I gas of 
SBS 0335-052 (with log (O/H); of the ionized gas equal to -4.70) and the ones seen in the H 1 
gas of other BCDs for which similar FUSE studies have been carried out, I Zw 18 [log (O/H); 
= -4.82], I Zw 36 [log (O/H); = -4.23], Mrk59 [log (O/H); = -4.01] and NGC 1705 [log (O/H); 
= -3.79]. In all five BCDs, log O/H of the neutral gas is in the range -5.0 - -4.5 (or [O/H] = 
—1.7 - -1.2), although the ionized gas of NGC 1705 is some 10 times more metal-rich than that 
of SBS 0335-052 and I Zw 18. This suggests that the metallicities of the neutral and ionized gas 
components are unrelated. We may speculate that the neutral gas in these five BCDs has been 
enriched to a common metallicity level of log (O/H) ~ -5 by population III stars. This speculation 
appears to be supported by observations of Lya absorbers which give similar oxygen abundances 
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in the intergalactic medium. 
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Table 1. Heavy element column densities (cm 2 ) in SBS 0335-052 



Species 


logiV (cm 2 ) a 


CLOUDY b 


H I (FUSE) 


2i.86tg:8i 


21.69lg-n 


C II 


17-5918:1 


14 KK+ - 02 
14.OO_ 002 


N I 


14.84+°f 


lO 1O+0.01 

±0. 1<J_Q QQ 


N II 


15 35 +i: ° 
io.oo_ 75 


I O yq + 0.02 

±0. IO_q Q1 


N II* 


13.70l° 3 25 70 




N in 


18.70±{ 8 ° 7 


15.29^8:18 


i 


i6.82i8:jj 


14.8618:81 


Si ii 


14 qn+ - 45 


io OK+0.04 

io.»o_ 004 


Ar I 


13-85l8| 5 


1 9 OO+0.01 

±Z.OZ_q Q1 


P ii 


14.0011:°, 


12.0618:81 


P ii* 


13.09l° 3 5 




Fe ii 


15 44 +ai 


13.6018:8? 



a Fitting with a single Voigt profile. Error 
bars are 1 a limits. 

b Total H I + H ii column density in the H 
ii region. 
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Table 2. Abundances in SBS 0335-052 used as input to CLOUDY 



Species 


log (iV(species)/A(H)) 


He 


-1.0915 a 


C 


-5.50 b 


N 


-6.26 a 





-4.67 a 


Ne 


-5.48 a 


Si 


-6.27 b 


P 


-8.05 c 


S 


-6.23 a 


Ar 


-6.93 a 


Fe 


-6.10 a 



a From Izotov et al. (1997). 

b From Izotov & Thuan (1999). 

c Adopting solar log P/O = -3.42 
from Grevesse & Noels (1996). 
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Table 3. Comparison between the observed and CLOUDY predicted optical line intensities 

normalized to H/3 



Ion 




Observed 3, 


CLOUDY 


3727 


[O ii] 


0.233 


n 91 7+0.046 
U - Z1 ' -0.041 


3835 


H9 


0.079 


u.u i »_q goo 


3868 


[Ne in] 


0.239 


U.ZOO_ 014 


3889 


He I + H8 


0.172 


n 104+0.002 


3968 


[Ne in] + H7 


0.239 


n 2 on+0.004 

u.zoy_ 003 


4101 


m 


0.255 


n 264+ - 000 
U.ZO4_ 000 


4340 


H 7 


0.476 


474+0-001 
u -^' ^-0.000 


4363 


[O in] 


0.109 


n ln2 +0.008 
U.1UZ_ 009 


4471 


He I 


0.035 


n35 +0 - 001 
u.uoo_ 000 


4658 


[Fe in] 


0.003 


02+ - 001 
u.uuz_ 000 


4740 


[Ar iv] 


0.009 


u - uu ' -0.001 


4861 


H/3 


1.000 


1 ooo +0 - 000 

i.uuu_ .ooo 


4959 


[O in] 


1.054 


1 nqq+0.061 
i.uyy_ 060 


5007 


[O in] 


3.155 


1 7 o+0.172 
°- 1 '°-0.175 


5876 


He i 


0.100 


n nq^+O.OOO 
u.uyo_ 000 


6300 


[O i] 


0.007 


02+ - 000 
u.uuz_ 001 


6312 


[S in] 


0.006 


n ni 0+0.000 

U.U10_oooi 


6563 


Ha 


2.745 


2 -° 4U -0.008 


6584 


[N ii] 


0.007 


u.uu^„o.ooo 


6678 


He 1 


0.027 


n27+ - 000 

v.vz, 1 _ .000 


6717 


[S ii] 


0.019 


n O+0.002 
u.u±o_ 003 


6731 


[S ii] 


0.017 


009+ - 002 
u.uuy_ 001 


7065 


He i 


0.039 


n n24+0- 001 

U.UZ4t_ 001 


7136 


[Ar in] 


0.014 


n29 +0 - 001 

u.u^y_ 001 



a From Izotov et al. (1997). 
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Table 4. CLOUDY predicted column densities in the H n region 



Ion log x a log N h 



C ii 


-1 64n +0 - 091 


u cc:+0.02 
14.OO_ 02 


N i 


9 OQE+0.108 


13 13+ - 01 

lo.lO_ 00 


N ii 


-1 7nn+ - 096 

L,uu -0 .102 


1 o 7O+0.02 
lO./O_ 01 


N in 


-0 141+o.on 

u - i4i -0.014 


i5-29±8:ig 


i 


-2 164+ - 100 


14-86±g;gi 


Si ii 


1 /Iko+0.073 


i3.96ig;8S 


P ii 


l.O»4_ 1Q9 


i2.06i8;81 


Ar I 


-2 436+ - 104 
z - 4OD -0.110 


1 9 oo+O.Ol 
iz -°^-0 .01 


Fe ii 


_1 qqi +0.098 

l.»»l_Q 1Q3 


13.6018:8? 



a Radially averaged ratio of the ion number to the total 
number of the element. 

b Column density in cnf 2 . 
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Fig. 1. — FUSE spectrum in the rest frame of SBS 0335-052 (shifted by 4050 km s _1 ). Here 
the data have been rebinned to a resolution of 0.2 A. Prominent interstellar absorption lines are 
indicated. The lines arising in SBS 0335-052 are marked on top and those in the Milky Way at 
bottom. In addition to the Hi Lyman series in SBS 0335-052, there are also strong interstellar 
absorption lines from C II, N I, N II, O I, Si n and Fe II. 
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Fig. 2. — Co-added profile of the Fen absorption line in SBS 0335-052. The profile and adjacent 
continua have been obtained by an error-weighted co-addition of 9 different Fe n lines at wavelengths 
A = 1055 A, 1063 A, 1064 A, 1082 A, 1097 A, 1122 A, 1125 A, 1143 A and 1145 A. The profile fits to 
individual lines have also been co-added, and the result is shown by a thick line. The dashed line 
shows the addition of the corresponding continua. Other lines can also be seen such as a H2 (J=l) 
line and a Fell line from the Milky Way at ~110 km s -1 . 
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Fig. 3. — Profile of the O I A1039A absorption line in the LiFla channel. The fit is shown by small 
dashes while the continuum is shown by long dashes. The other absorption line at A ~ 1053.4 A is 
due to Galactic H2. Its fit is also shown by small dashes. The sum of the two profiles is shown by 
the thick solid line which fits well the data (histogram). 
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- O Column densities in the H II region 

Element 

Fig. 4. — Heavy element column densities in the H I and H n gas of SBS 0335-052. The filled 
circles show the column densities derived from the FUSE data. The open circles show the column 
densities in the Hn region as calculated from CLOUDY modeling of the optical emission lines. 
The column densities in the H I cloud follow the same pattern as the ones in the H n gas, although 
they are larger by 1-2 orders of magnitude. While this pattern reflects the cosmic abundances of 
the considered species, the fact that we recover this pattern shows the consistency of the measured 
column densities in the Hi gas although they have been derived from low S/N data. The lower 
column densities in the ionized gas show that the absorption lines in the FUSE spectra arise in the 
neutral interstellar gas. 
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Fig. 5. — Relative ratios of FUSE column densities (filled circles) in the H I gas compared to the 
corresponding abundance ratios in the H n gas (stars) in SBS 0335-052. For each FUSE column 
density, the ionization state shown is the main one in the H I region and the ionization correction 
factors are very near unity. For the H n region, total element abundances including all ionization 
states are plotted. From left to right are the C/O, N/O, Si/O, Ar/O, Fe/O and O/H ratios. The 
right ordinate is for the O/H ratio, while the left ordinate is for all other abundance ratios. Note 
in particular that the O/H ratios in the Hi and Hll regions are very similar. 
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